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We present the results of detailed theoretical investigations of changes in local density of total 
electronic surface states in 2D anisotropic atomic semiconductor lattice in vicinity of impurity atom 
for a wide range of applied bias voltage. We have found that taking into account changes in density 
of continuous spectrum states leads to the formation of a downfall at the particular value of applied 
voltage when we are interested in the density of states above the impurity atom or even to a series 
of downfalls for the fixed value of the distance from the impurity. The behaviour of local density of 
states with increasing of the distance from impurity along the chain differs from behaviour in the 
direction perpendicular to the chain. 
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I. INTRODUCTION 

Influence of different impurities on the semiconduc- 
tor local density of surface states was widely studied 
experimentally and theoretically. Most of the experi- 
ments were carried out with the help of scanning tun- 
neling microscopy/spectroscopy technique [l|, (2|, Q. 
Theoretical investigations of single impurities and clus- 
ters influence on density of surface states deals with 
Green's functions formalism or based on the 

total-energy density-functional calculations using first- 
principle pseudo-potential [g]. Numerical calculations 
based on the tight-binding model are also carried out 
0- 

Most of the theoretical calculations don't take into ac- 
count modification of local density of continuous spec- 
trum states due to the influence of impurity atom on 
semiconductor surface which we consider to play an im- 
portant role in the formation of peculiarities in the local 
density of total surface states . So in the present work 
we suggest a simple model of anisotropic atomic lattice 
with impurity atom and we pay special attention to the 
changes of local density of continuous spectrum states. 
This model suits well for theoretical investigation of 7T- 
bonded chains on the reconstructed Ge or Si surfaces 
Q. It can be also used for investigation of the sublat- 
tices on the cleaved planes of Am By semiconductors. 
We have found that the view of local density of surface 
states (amount of downfalls and their shape) strongly 
differs depending on the value of the distance from the 
impurity atom position and from the direction of obser- 
vation (along the atomic chain or perpendicular to the 
atomic chain). It will be shown that downfall in the res- 



onance when we are interested in the density of states 
above the impurity atom can transforms to a series of 
downfalls or even to a peak for different values of the 
distance from the impurity. 



II. THE SUGGESTED MODEL AND MAIN 
RESULTS 

We shall analyze 2D anisotropic atomic lattice formed 
by the similar atoms with energy levels e\ and similar 
tunneling transfer amplitudes between the atoms t along 
the atomic chain. The interaction between atomic chains 
is described by tunneling amplitude T, which has the 
same value for all the similar atoms in the chain (Fig. 1). 
Distance between the atoms in the atomic chain is equal 
to a, distance between the atoms in the neighboring 
chains is equal to b. Atomic lattice includes impurity 
atom with energy level £<j, tunneling transfer amplitude 
from impurity atom to the nearest atoms in the atomic 
chain r and to the nearest atoms in the neighbor chains 
3 . 

The model system can be described by the Hamilto- 
nian: H: 



H — H + H imp + H t 

Ho = ^ Sicfci + ^ tc t C 3 + X] T °k Cl + h- C - 
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H imp 
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Ho is a typical Hamiltonian for atomic lattice with hop- 
pings without any impurities. H tU n describes transitions 
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between impurity atom and neighboring atoms of the 
atomic lattice. Hi mp corresponds to the electrons in the 
localized state formed by the impurity atom in the atomic 
chain. 

Indexes i,j correspond to the direction along the chain; 
indexes k,l correspond to the direction perpendicular to 
the chain. 
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FIG. 1: Schematic diagram of 2D atomic lattice with impu- 
rity atom. 

We shall use diagramm technique in our investigation 
of atomic chain local density of states. 

The dependence of local density of states on the dis- 
tance along the atomic chain and in the direction perpen- 
dicular to the atomic chain in the presence of impurity 
atom is described by the equation: 

p(u, f) = —Sp( Im J2 G R {^ «1, u)e iHf! e iKV ) (2) 



Where f = (x,y), k = (k x ,k v ) and K\ = (ri xl ,n yl ). 
Green function G r (k, k[,u>) corresponds to the electron 
transition from the impurity to the semiconductor con- 
tinuum states and can be found from the system of equa- 
tions: 



u dd 



rtR _ riOR ^.r-iR , s-fOR. rp \ "* r R. 
< - r K x 0drf — '~ r K, m OK : cO T{ - r dd T ^OKjO- 1 / , ^k.k, 

riR 

^OKydd 

r R _ r oR , r an r R , r QRc* ST r R 

u dd — u dd ' u dd T 2_^i ^KxOdd ' u dd ^ ^Okj, 



dd 



dd 



U Kdd — L 'KK TLT K :c 0dd "+" U KK JLr Of 



K v dd 



n R _ s-tOR , rtOR^rtR , rtVR^r.R, 

U KK! — u itFZ "+~ PIk ddK± T ^KK^WciKl 



Where zero Green function is evaluated for the 2D 
atomic lattice without any impurities and has the form: 

qOR / \ _ \ 

KlK » l ; oj - £i - 2t ■ cos(k x a) - 2T • cos(k y b) 

(3) 

Substituting the expression for Green function 
G R (n, Ki, u>) obtained from system into equation (|2|) and 
performing summarization over wave vectors k x and k x \ 



(k y and k y \) we get the final expression for the local den- 
sity of continuous spectrum states along (perpendicular) 
the atomic chain p vo i um e(x) (p V oiume(y))- 



Pvolume{u,x) = Po(uj) 



{LQ-e d f + 7 2 •(!- f(2k x (u)x) 
(cj-e d ) 2 +7 2 



I , r \ ^-e d f+ 1 2 -{l~f{2k y {uj)y) lA ^ 

(w - £ d y + 7 

Where parameter 7 = (r 2 + 3 2 ) • /0o(w) corresponds 
to relaxation rate of electron distribution at the local- 
ized state formed by impurity atom, po(u>) is a local den- 
sity of states for the atomic chain without any impuri- 
ties. Functions f(2k x (uj)x) and f(2k y (uj)y) are periodi- 
cal and have the property: f(2k x (uj)x) = 1 if x = and 
f(2k y (oj)y) = 1 if y = 0. If both directions are equiv- 
alent f(2k x (uj)) = f(2k y (uj)y) = Jo {2k x {u))). Expression 
for k x (u) or k y (w) can be found from the dispersion law 
of the 2D atomic lattice which has the form. 



L)(k x ,k y ) — 2t ■ cos(k x a) + 2T ■ cos(k y b) 



(5) 



Impurity atom density of states has lorentzian form 
line shape and can be evaluated as: 



Pimpurity{u) — ■ IfU N ^ G R A(J) 

IT ' * 



{u-e d ) 2 + 1 2 



(0) 



Local density of total surface states is the result of sum- 
marization between local density of continuous spectrum 
states and impurity atom density of states. 



p(uj) — pvolume(^) + Pimpurity(oj) 



(7) 



Let's start from the ID case of the atomic chain. In this 
case it is necessary to put in the Hamiltonian: b = 0, 
T = and 3 = 0. Final expression for the local density 
of continuous spectrum states p V oiume will have the form: 



Pvolumeiu) — Po(ti>) 



(uj - e d ) 2 + 7 2 • (1 - cos(2fc(w)r) 
(w - e d ) 2 + 7 2 



(8) 



Expression for k(u>) can be found from the dispersion law 
of the ID atomic chain. 

Typical numerical results for local density of total 
surface states and local density of continuous spectrum 
states calculated above the impurity atom in ID case 
(distance value is equal to zero (r — 0)) are shown on 
(Fig. 2a, b). For the local density of continuous spectrum 
states (Fig. 2a) a downfall exist in the resonance when 
energy is equal to the impurity atom energy level depo- 
sition (ui = e d ). Width of the downfall depends on the 
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FIG. 2: a) Local density of continuous spectrum states in the case of the distance from the impurity atom along the atomic 
chain equal to zero, b) Local density of surface spectrum states in the case of the distance from the impurity atom along the 
atomic chain equal to zero. c)-f) Local density of continuous spectrum states (black line) and local density of surface spectrum 
states (grey line) for the different values of the distance from the impurity atom along the atomic chain. For all the figures 
values of the parameters a — 1, t — 1, 5, Ed = 0, 6 are the same. 



parameters of the atomic chain, such as relaxation rate 
or tunneling transfer amplitude, it rises with the increas- 
ing of tunneling transfer amplitude from impurity atom 
to the neighbor atoms of the atomic chain. Local den- 
sity of surface states (Fig. 2b) has lorentzian form with 
a downfall in the resonance. With the increasing of re- 
laxation rate (increasing of r) the downfall depth at the 
top of the peak decreases, resonance peak shape spreads 
and it's amplitude falls down. 

Now let's start to analyze the dependence of local den- 
sity of continuous spectrum states and local density of 
surface states at the fixed value of the distance r along the 
atomic chain from the impurity atom position (Fig. 2c-f). 
We shall again start from the local density of continuous 
spectrum states (black lines on Fig. 2). 

When the value of a distance is not equal to zero a se- 
ries of downfalls in the local density of continuous spec- 
trum exists. Amount of downfalls increases with the in- 
creasing of distance value and downfalls amplitude de- 
creases when energy aspire to the edges of the band. The 
most significant amplitude of the downfalls corresponds 
to the vicinity of the resonance region. It is clearly evi- 
dent that positions of the downfalls on the energy scale 
can be found from the equation 27m = 2k(u>)r where 
n is an integer number. This means that numerator of 



the equation ( [5]) is equal to zero. When the distance is 
not equal to zero not only a downfall in the resonance 
(Fig. 2d,f) but also a peak (Fig. 2c, e) can exist in the 
local density of continuous spectrum states. At the fixed 
parameters of the atomic chain existance of a downfall 
or a peak in the resonance is determined by the value 
of the distance. Local density of surface states is shown 
by the grey line on Fig. 2c-f. Comparison between local 
density of surface states and local density of continuous 
spectrum states for the fixed value of r make it clearly 
evident that impurity atom density of states can drasti- 
cally influence on the local density of continuous spec- 
trum states. Result depends on the value of tunneling 
transfer amplitude from impurity atom to the nearest 
atoms of the chain. We have found that the most signifi- 
cant influence corresponds to the situation when tunnel- 
ing transfer amplitude between the atoms in the chain 
t significantly exceeds tunneling transfer amplitude from 
impurity atom to the atoms of the chain r. In this case 
downfall in the resonance becomes a peak with a small 
downfall on the top of the peak (Fig. 2e). So only one sig- 
nificant downfall exists in the resonance region and there 
are no downfalls at the energies different from the reso- 
nance value (Fig. 2c)in the local density of surface states, 
while in the continuous spectrum density of states a series 
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FIG. 3: a) Local density of continuous spectrum states (black line) and local density of surface spectrum states (grey line) for 
the different values of the distance from the impurity atom along the atomic chain. For all the figures values of the parameters 
a = 1, b = 2, t = 1, 5, T = 1, 2, r = 0, 6, 9 = 0, 3, £d = 0, 6 are the same. 



of downfalls exists. With increasing of r local density of 
surface states differs from the local density of continuous 
spectrum states only by the amplitude of the downfalls 
(Fig. 2e,f). In this case number of downfalls is the same 
in comparison with the continuous spectrum density of 
states, downfalls don't change their shape or position on 
the energy scale and contribution to the local density of 
surface states from the impurity atom can be considered 
to be a background. This effect can be qualitatively un- 
derstood by the following way: with increasing of transfer 
amplitude from impurity atom to the atom of the chain 
increases relaxation rate and lorentzian form peak of im- 
purity atom density of states spreads and it's amplitude 
decreases. 

Now let's start to analyze 2D atomic lattice. Nu- 
merical results for local density of total surface states 
and local density of continuous spectrum states calcu- 
lated above the impurity atom in perpendicular direc- 
tions (along the atomic chain and perpendicular to the 
atomic chain) are shown on Fig. 3a, Fig. 4a. It's clear 
that in this case all the results are equal for both direc- 
tions and downfall in the continuous spectrum density 
of states or a peak in the local density of surface states 
poses just the same properties as in the case of ID atomic 
chain. 

Let's analyze the situation when the value of the dis- 



tances from the impurity atom position along the atomic 
chain (Fig. 3b-f) and perpendicular to the atomic chain 
(Fig. 4b-f) are not equal to zero. We shall start from the 
local density of continuous spectrum states (black lines 
on Fig. 3,4). 

In this case number of downfalls, their shape and po- 
sition on the energy scale in each of the perpendicular 
directions can be found from the dispersion law just in 
the same way as for II? atomic chain. When the distance 
is not equal to zero not only a downfall in the resonance 
(Fig. 3d,e; Fig. 4e,f) but also a peak (Fig. 3c, f; Fig. 4c, d) 
can exists in the local density of continuous spectrum 
states. 

We have found distance interval for both directions 
(along the atomic chain and perpendicular to the atomic 
chain) where exists replacement of a peak by a down- 
fall (Fig. 3c-f; Fig. 4c-f) in local density of continuous 
spectrum states (and also in the local density of sur- 
face states). Moreover peak in one direction can cor- 
responds to a downfall in the another direction. Let's 
analyze this interval carefully. We shall start from the 
distance value when peaks in the resonance in local den- 
sity of continuous spectrum states for both directions can 
be seen (x(y) — 15 -a) (Fig. 3c; Fig. 4c). With increasing 
of the distance value (x(y) = 20 • a) a resonance peak in 
the direction perpendicular to the atomic chain still ex- 
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FIG. 4: Local density of continuous spectrum states (black line) and local density of surface spectrum states (grey line) for 
the different values of the distance from the impurity atom perpendicular to the atomic chain. For all the figures values of the 
parameters a = 1, b = 2, t = 1, 5, T = 1, 2, r = 0, 6, = 0, 3, Sd = 0, 6 are the same. 



ists (Fig. 4d) and in the direction along the atomic chain 
a downfall appears (Fig. 3d). Further increasing of the 
distance value (x(y) = 25 • a) shows that in the direction 
along the atomic chain a downfall still exists (Fig. 3e) 
and in the perpendicular direction a downfall substitutes 
peak (Fig. 4e). Finally when the value of distances in 
both directions becomes equal to x(y) = 30 • a in the per- 
pendicular direction a downfall still exists (Fig. 4f) and 
in the direction along the atomic chain a resonance peak 
can be seen (Fig. 3f). 

Local density of total surface states is shown by the 
grey lines on Fig. 3,4. It is clearly evident that for the 
studied parameters of the system taking into account im- 
purity atom density of states slightly changes local den- 
sity of total surface states in comparison with local den- 
sity of continuous spectrum states for both directions. In 
this case number of downfalls is the same in comparison 
with the continuous spectrum density of states, downfalls 
don't change their shape or position on the energy scale. 

III. CONCLUSION 

In this work we have shown that taking into account 
changes of the local density of continuous spectrum states 



formed by the presence of the impurity atom in the 2D 
anisotropic atomic lattice or even in the ID atomic chain 
leads to significant modification of the total local density 
of surface states and consequently to the modification of 
STS spectra. We have found that a downfall exists in 
the STS spectra measured just above the impurity when 
impurity atom energy level is equal to the applied bias 
voltage. With changing of the distance from the impurity 
a series of downfalls is formed on the energy scale both in 
the local density of continuous spectrum states and in the 
local density of total surface states. Number of downfalls 
and their position are determined by the atomic lattice 
dispersion law. It was shown that at some values of the 
distance from the impurity a peak can exist in the reso- 
nance region instead of a downfall. We have found that 
behaivour of local density of surface states depends on 
the direction of the observation. Switching on and off 
of impurity atom in both directions was found. This ef- 
fect can be well observed experimentally with the help of 
STM/STS technique. 
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